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C and N Abundances in Stars At the Base of the Red Giant Branch in M5 ^ 


Judith G. Cohen^, Michael M. Briley^ and Peter B. Stetson^’^’® 

ABSTRACT 

We present an analysis of a large sample of moderate resolntion Keck LRIS spectra 
of snbgiant {V ~ 17.2) and fainter stars in the Galactic globular cluster M5 (NGC 5904) 
with the goal of deriving C and N abundances. Star-to-star stochastic variations with 
significant range in both [G/Fe] and [N/Fe] are found at all luminosities extending to 
the bottom of the RGB at My ~ +3. Similar variations in GH appear to be present in 
the main sequence turnoff spectra, but the signal in the cnrrent sample is too low for 
a detailed analysis. The variations seen among the M5 subgiants are consistent with 
the abundances found earlier by Briley et al. (1992) for brighter giants in this cluster. 

There is thus no sign of a change in the behavior of G and N with evolutionary stage 
over the full range in luminosity of the RGB and SGB, although a systematic decrease 
with luminosity in the mean [G/H] smaller than a factor of 2 cannot be ruled out with 
confidence at present. 

The C and N abundances appear strongly anti-correlated, as would be expected 
from the GN-cycle processing of stellar material. Yet the present stars are considerably 
fainter than the RGB bump, the point at which deep mixing is believed to set in. 

On this basis, while the observed abundance pattern is consistent with proton capture 
nucleosynthesis, we infer that the site of the reactions is likely not within the present 
sample, but rather in a population of more massive (2 - 5Mq) now defunct stars. 

The range of variation of the N abundances is very large and the sum of C-l-N 
increases as C decreases. To reproduce this requires the incorporation not only of 
CN but also of ON-processed material. Furthermore, the existence of this correlation 
is quite difficult to reproduce with an external mechanism such as “pollution” with 
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material processed in a more massive AGB star, which mechanism is fundamentally 
stochastic in nature. We therefore suggest that although the internal mixing hypothesis 
has serious flaws, new theoretical insights are needed and it should not be ruled out yet. 


Subject headings: globular clusters: general — globular clusters: individual (M5) — 
stars: evolution - stars: abundances 

1. Introduction 

By virtue of their large populations of coeval stars, the Galactic globular clusters present us 
with a unique laboratory for the study of the evolution of low mass stars. The combination of their 
extreme ages, compositions and dynamics also allows us a glimpse at the early history of the Milky 
Way and the processes operating during its formation. These aspects become even more significant 
in the context of the star-to-star light element inhomogeneities found among red giants in every 
cluster studied to date. The large differences in the surface abundances of G, N, O, and often Na, 
Mg, and A1 have defied a comprehensive explanation in the three decades since their discovery. 

Proposed origins of the inhomogeneities typically break down into two scenarios: 1) As G, N, 
O, Na, Mg, and A1 are related to proton capture processes at CN and CNO-burning temperatures, 
material cycled through a region in the upper layers of the H-burning shell in evolving cluster 
giants may be brought to the surface with accompanying changes in composition. While standard 
models of low mass stars do not predict this “deep mixing”, several theoretical mechanisms have 
been proposed (e.g., the meridional mixing of Sweigart & Mengel 1979, and turbulent diffusion, 
Charbonnel 1994, 1995) with varying degrees of success. Moreover, there is ample observational 
evidence that deep mixing does take place during the red giant branch (RGB) ascent of metal-poor 
cluster stars (see the reviews of Kraft 1994 and Pinsonneault 1997 and references therein). 2) It 
has also become apparent that at least some component of these abundance variations must be in 
place before some cluster stars reach the giant branch. Spectroscopic observations of main sequence 
turn-off stars in 47 Tuc (Hesser 1978; Hesser &: Bell 1980; Bell, Hesser, & Cannon 1983; Briley, 
Hesser, &: Bell 1991; Briley et a/.1994, 1996; Cannon et al. 1998) and NGC 6752 (Suntzeff & Smith 
1991, Gratton et al. 2000), as well as our own work in M71 (Cohen 1999, Briley & Cohen 2001, 
Ramirez &: Cohen 2002) have shown variations in CN and CH-band and Na line strengths consistent 
with patterns found among the evolved giants of these clusters. The assumption that these low 
mass cluster stars are incapable of both deep dredge-up and significant CNO nucleosynthesis while 
on the main sequence leads to the possibility that the early cluster material was at least partially 
inhomogeneous in these elements or that some form of modification of these elements has taken 
place within the cluster. Suggested culprits include mass-loss from intermediate mass asymptotic 
giant branch stars and supernovae ejecta (see Cannon et al. 1998 for an excellent discussion of these 
possibilities). 
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Thus the observed light element inhomogeneities imply that there is some aspect of the struc¬ 
ture of the evolving cluster giants which remains poorly understood (the deep mixing mechanism), 
that the early proto-clusters may have been far less homogeneous, that intermediate mass stars 
may have played a greater role in setting the composition of the present day low mass stars than 
previously thought, etc. Indeed, the evidence cited in the reviews above has led many investigators 
to suggest that a combination of processes is responsible, i.e., many clusters contain star-to-star 
inhomogeneities established early in their histories which have subsequently been further altered 
by deep mixing during the ascent of the RGB. This of course greatly exacerbates the difficulty of 
achieving an understanding of these issues, as a knowledge of the composition of the more easily 
observed bright red giants will not tell the whole story of their chemical history - one must also 
understand the makeup of the main sequence stars. 

In the present paper, we continue our earlier work on M71 by exploring the CH and CN 
band strengths in a sample of low luminosity stars in the somewhat more metal poor globular 
cluster M5. We adopt current values from the on-line database of Harris (1996) for the apparent 
distance modulus of M5 at V of 14.31 mag with a reddening of E(B-V) = 0.03 mag. Recent CCD 
photometric studies of this cluster, focusing primarily on its age, are given by Johnson & Bolte 
(1998) and Stetson et al. (1999). Sandquist et al. (1996) discuss the predominantly blue horizontal 
branch of M5. We adopt the metallicity [Fe/H] = —1.21 dex found by Ivans et al. (2001) in a high 
dispersion abundance analysis of a large sample of stars on the upper giant branch of M5. 

We describe the sample in ^ and ^ We outline our measurement of the molecular band indices 
and their interpretation in terms of the scatter in §|^. With an assumption about the O abundance, 
these are converted into C and N abundances, from which we find a strong anti-correlation between 
C and N in §1^. A discussion of our results together with a comparison with the trends seen among 
the red giants in M5 and in other globular clusters is given in and A brief summary concludes 
the paper. 


2. Photometric Database 

The photometry of M5 employed here was carried out as part of a larger program to provide 
homogeneous photometry for star clusters and nearby resolved galaxies (Stetson, Hesser & Smecker- 
Hane 1998, Stetson 2000). In the present instance, the photometry is based on data obtained 
during 18 observing runs between 1984 and 1998. Of these, some observing runs were carried out 
by PBS, data from other runs were donated by collaborators (M. Bolte, J. Hesser, R. McClure, 
N. Suntzeff), and data from still other runs were obtained through the services of the Isaac Newton 
Group Archive. Data origins include the 4m and 0.9m telescopes on Kitt Peak and Cerro Tololo, the 
Canada-Prance-Hawaii Telescope, the Isaac Newton Telescope, and the Jacobus Kapteyn Telescope. 
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[[] In all, some 374 images of parts of M5 in filters corresponding to the BVI bandpasses were 
analyzed. Since different field sizes and different field centers were available from the various 
observing runs, no individual star actually appeared in more than 226 of those images. 

Employing only data obtained under photometric conditions on nights when numerous obser¬ 
vations of fundamental and secondary standards (see Stetson 2000) were also obtained, we defined 
a network of 649 local photometric standards in the M5 field. A star was considered a “local 
standard” in a given photometric bandpass only if it had been measured in that bandpass under 
photometric conditions on at least five occasions, if the standard error of the mean magnitude in 
that bandpass was less than 0.02 mag, and if the star showed no evidence for intrinsic variability 
greater than 0.05 mag, root-mean-square, in excess of random measuring errors when all bandpasses 
were considered together. Once this network of local standards had been created, it was possible to 
use them to calibrate the zero points of frames taken under non-photometric conditions (the color 
terms of the photometric transformations having been obtained through observations of fundamen¬ 
tal standard fields), allowing us to include these latter frames in the overall photometric solution. 
Adding images taken under non-photometric conditions does nothing to improve the accuracy of 
our photometric transformations to the fundamental standard system in absolute terms; however, 
it does allow us to increase the precision of the photometry of individual stars relative to the mean 
system defined by our 649 local photometric standards. 

In our experience, photometry from datasets such as those employed here typically display 
an external accuracy of roughly 0.02 mag per observation; this level of observation-to-observation 
scatter is probably dominated by temporal and spatial fluctuations in the instantaneous atmospheric 
extinction, and probably also by the difficulty of obtaining truly appropriate flat-field corrections in 
the presence of such effects as scattered light, ghosts, fringing and spectral mismatch between the 
flat-held illumination and the astronomical scene. In the present instance, our absolute photometry 
for the 649 local standards is based upon a median number of 10, 29, and 19 observations in the 
B, V, and I bandpasses. Taking 0.02 mag per observation as a typical external uncertainty, we 
therefore expect that our photometry is referred to the fundamental system (i.e., that of Landolt 
1992) with an absolute accuracy of order 0.006, 0.004, and 0.005 mag, standard error of the mean, 
in B, V, and I. Relative magnitude and color differences among different stars in our sample can be 
smaller than this, depending upon their brightness and crowding and the number of observations 
that went into each determination. 

The absolute astrometry of our catalog is based upon the United States Naval Observatory 
Guide Star Catalogue I (A V2.0; henceforth USNOGSC, Monet et al. 1998), access to which is 
obtained by PBS through the services of the Canadian Astronomy Data Centre. The authors of 
the USNOGSC claim a typical astrometric error of 0.25arcsec, which they believe is dominated 


^Although some Keck images are available, very short exposure times of only 1 or 2 seconds were used, and the 
photometric accuracy of the data is then compromised at the level of 2% due to issues involving the shutter open/close 
time and its uniformity over the field. In addition, the field of the available Keck images is small. 
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by systematic errors in the calibration procedure. When transforming relative (x, y) positions 
from large-format CCD images to absolute right ascensions and declinations from the USNOGSC 
for stars in common, we typically find root-mean-square differences of 0.3 to 0.4arcsec in each 
coordinate. Some of this is clearly due to proper-motion displacements accumulated during the 
forty-plus years between the obtaining of the first Palomar Observatory Sky Survey and our CCD 
data. However, a significant part of the differences is also due to the lower angular resolution of 
the Schmidt plates as scanned by the Precision Measuring Machine (built and operated by the 
U. S. Naval Observatory Flagstaff): particularly in crowded helds such as the outskirts of globular 
clusters, a single entry in the USNOGSC is occasionally found to correspond to the photocenter of 
a close pair or a clump of several stars in the CCD imagery. 

As a result of these non-Gaussian errors (i.e., proper motions and blending), we perform our 
astrometric transformations using an iterative procedure wherein 20-constant cubic polynomials 
are used to approximate the transformation of the (x, y) positions measured in our GGD images to 
standard coordinates obtained from a gnomonic projection of the right ascensions and declinations 
listed in the USNOGSC. After each iterative step, stars lying within a certain radial distance of 
their predicted positions are used to obtain an improved geometric transformation with individual 
weights gradually tapering from 1.0 for perfect positional agreement to 0.0 for stars lying exactly 
one critical radius from their predicted positions; stars lying farther than one critical radius from 
their predicted positions are ignored altogether. The critical radius distinguishing acceptable cross¬ 
identifications starts out fairly generous (~ 5 arcsec), but then this radius is gradually reduced until 
it is down to a value of one arcsecond; the iterative process then continues, maintaining this critical 
radius of one arcsecond, until a stable list of cross-identifications and transformation constants 
is achieved. As an indication of the percentage of entries in the USNOGSG whose positions we 
consider no longer strictly valid, due either to proper motions or to blending, we find that 4,709 
USNOGSC entries lie within 5 arcsec of one of our CCD detections; 4,652 lie within 2 arcsec of a 
CCD detection; and 3,886 he within 1 arcsec. When the sample is restricted to these 3,886 cross¬ 
identifications, the unweighted root-mean-square positional difference is approximately 0.34 arcsec 
in both X and y. Thus, what is essentially a 3 -cj clip results in a rejection of some 17% of the 
“possible” cross-identifications that agree to within 5 arcsec. 

Of these 3,886 acceptable cross-identifications, only a dozen or so lie within 200 arcsec of the 
cluster center; none at all lie within 160 arcsec. Therefore it is possible that interpolation of the 
astrometric transformations across this gap could be subject to small systematic errors, probably 
not as large as 0.1 arcsec. Throughout the region of our field that is well populated by USNOGSG 
stars (including essentially all of the stars in our present spectroscopic sample), we expect systematic 
errors of our right ascensions and declinations on the system of the USNOGSC to be of order 
0.34 arcsec/\/3886 — 20, or ~ 0.01 arcsec. We have no independent information on the accuracy 
with which the USNOGSC coordinate system is referred to an absolute inertial reference frame 
of right ascension and declination. Individual random errors in our coordinate measurements are 
probably not much better than 0.02 arcsec on a star-by-star basis, the errors becoming somewhat 
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worse than this for the fainter and more crowded stars in our photometric/astrometric sample. 

Stars are identified in this paper by a name derived from their J2000 coordinates, so that star 
C12345_5432 has coordinates 15 12 34.5 +2 54 32. 


3. Spectroscopic Observations 

The initial sample of stars consisted of those from the photometric database located more 
than 150 arcsec from the center of M5 (to avoid crowding) with 16.9 < V < 17.35 and with 
0.86 < (y — I) < 0.96, i.e. subgiants at the base of the RGB. (A preliminary version of the 
photometric catalog described in was used for this purpose.) From this list, two slitmasks 
containing about 25 slitlets each were designed using JGC’s software. The center of the first field 
was roughly 3.2 arcmin E and 0.5 arcmin S of the center of M5, while the center of the second field 
was located roughly 2.5 arcmin W and 0.5 arcmin S of the cluster center. 

These slitmasks were used with LRIS ( pke et al. 1995| ) at the Keck Observatory in May 2001; 
three 800 sec exposures were obtained with each slitmask. The exposures were dithered by moving 
the stars along the length of the slitlets by 2 arcsec between each exposure. Because of the crowded 
fields, there were often more than one suitably bright object in each slitlet. Hence subtraction of 
sequential exposures was not possible, and they were reduced individually using Figaro (Shortridge 
1988) scripts, then the ID spectra for each object were summed. The width of the slitlets was 
0.7 arcsec, narrower than normal to enhance the spectral resolution. LRIS-B was used with a 600 
line grism giving a dispersion of 1.0 A/pixel (3.0 A resolution for a 0.7 arcsec wide slit). This gave 
good coverage of the region from 3600 to 4800 A, including the key GN band at 3885 A and the 
G band of GH at 4300 A. Technical problems related to the layout of the slitmasks led to some 
of the spectra not extending blue enough to reach the uv GN band. The red side of LRIS was 
configured to use a 1200 g/mm grating centered at Ha with the intention of providing higher 
accuracy radial velocities. The dispersion is then 0.64 A/pixel (29 km s“^/pixel) or 1.9 A/spectral 
resolution element. Unfortunately, the LRIS-R shutter was not functioning on the first night, and 
so only one of the two slitmasks used has matching Ha spectroscopy. 

In addition to the primary sample described above, since these fields are rather crowded, other 
stars sometimes serendipitously fell into the slitlets, and their spectra were also reduced. We refer 
to the latter as the secondary sample. As might be expected from the luminosity function, most 
of the secondary sample consists of stars at or just below the main sequence turnoff. Two stars in 
the secondary sample (C18386_0713 and C18465_0214) are early type stars. The former is almost 
certainly a BHB in M5, while the status of the latter is unclear; we subsequently ignore both of 
them. 
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3.1. Membership 

Given the high galactic latitude of M5 (b = 46.8°) and the proximity of our fields to the center 
of the cluster, we expect minimal field star contamination among the primary sample of stars. 
However, the secondary sample selection is based solely on spatial position and hence may contain 
a higher fraction of field stars. There are three indicators we use to establish membership. The 
first is deviation from the cluster sequences in a color-magnitude diagram. Stars which are off the 
M5 color locus in the {V,I) CMD by more than 0.1 mag were considered possible non-members. 
The second is the strength of the absorption features in the LRIS multi-slit spectra, ignoring the 
molecular bands. The metallicity of M5 is sufficiently low that near solar metallicity field stars of 
similar colors are easily distinguished from cluster members. We also examine the radial velocity 
by cross correlating the red spectra over the regime from 6400 to 6620 A. 

The radial velocity of M5 from the compilation of Pryor &: Meylan (1993) is -1-53.1 km s“^, with 
a = 4.9 km s“^. The histogram of heliocentric radial velocities we have measured for 37 stars is 
given in Figure [^. Only one star has a radial velocity inconsistent with membership. Rejecting this 
very discrepant object, with 36 remaining stars, we find a mean Vr = -1-52.8 km s“^, with a = 11.8 
km s“^. This value of a corresponds, after removing the intrinsic stellar velocity dispersion, to an 
instrumental error of 0.4 pixels in the focal plane CCD detector of LRIS-R, which seems reasonable. 

Table gives a listing of those stars in our sample which are probably not members based on 
these three tests. The final column gives the distance in {V — I) color at the fixed V mag of the star 
between the color of the star and that of the cluster locus. There are two dehnite non-members 
which appear to be more metal rich field stars. There are three stars which are off the cluster (V, I) 
color locus by slightly more than the adopted cutoff of 0.1 mag, but which we believe to be cluster 
members; all of them are crowded with a star of comparable brightness within 2 arcsec, and all are 
from the primary sample. There is also one faint star from the secondary sample which is either a 
blue straggler in M5 or an early type background object, and which we subsequently ignore. 

Spectra taken with the blue-side of LRIS are available for all the stars, but at a lower dispersion. 
After excluding the probable non-members listed in Table ||, cross correlations over the range 3800 
- 4400 A show the same la rms dispersion as does the red side, i.e. 0.4 pixels, here equivalent to 
29 km s“^. While the accuracy of these radial velocities, given the low spectral resolution, is not 
sufficient to confirm membership for any particular star, the small dispersion does demonstrate 
that almost all of the sample studied here must be members of M5. 


3.2. The Sample of Probable Members of M5 

There are 59 stars left in our M5 sample after the probable non-members and early type stars 
have been eliminated. In Figure ^ these are shown superimposed on the (F, I) CMD diagram, 
where a decreasing fraction of stars in the photometric database are plotted as V becomes fainter 


to maintain clarity in the figure. All but one of these are at least 1.5 mag fainter in V than the HB 
in M5; there is no possible confusion with AGB stars at these low luminosities. 

Figure ^ shows the region of the 4300 A G band of CH in the spectra of two of the stars in the 
primary sample. These stars have essentially the same stellar parameters {T^ff and log(( 7 )) lying 
at about the same place in the cluster CMD, yet their G bands differ strongly. The GN band near 
4200 A is too weak to be used in these metal poor stars; the ultraviolet GN band near 3885 A must 
be used instead. From this figure alone, we can anticipate one of the key results of our work, the 
large scatter in C abundance we will find among M5 members at the base of the RGB at V ~ 17.2 
mag. We can also see that the GN bands are much weaker, and a careful measurement will be 
required, which is complicated by the difficulty of measuring anything resembling a continuum at 
or near 3885 A. 

Table ^ lists the photometry for the sample of 59 stars (plus one BHB star) which we believe 
are members of M5. 


4. Measurement of CH and CN Indices 

For each spectrum, S(3839) and I(GH) indices sensitive to absorption by the 3885 A CN band 
and the 4300 A CH band respectively (see Briley &: Cohen 2001) were measured. These indices are 
listed in Table ^ and plotted in Figure |^. The error bars (drawn at 2 ct) have been calculated strictly 
from Poisson statistics based on the signal present in the feature and continuum bandpasses. 

It is extremely difficult to flux spectra taken through slitmasks because of the varying slit losses 
and the possibility of atmospheric dispersion affecting the spectra. The latter issue is of somewhat 
more concern than usual to us since we are working in the blue and UV with a narrow (0.7 arcsec 
wide) slit. Carrying out the observations with the length of the slit set to the parallactic angle, 
which is the usual method for eliminating atmospheric dispersion for single slit observations, cannot 
be used for multislit observations as the position angle is fixed by the design of the slitmask. It is 
for these reasons that no attempt was made to flux the spectra. 

The reduction process for the G band index of CH, whose feature bandpass is 4285 - 4315 A (ad¬ 
justed by 0.76 A for the average radial velocities of the present stars), is independent of this suite 
of issues as a continuum can be established both on the blue and on the red side of the feature. 
However, it is extremely difficult to establish a continuum on the blue side of the 3885 A CN band, 
and hence we are forced to rely on a single sided index, with a continuum bandpass determined 
only redward of the feature of interest. 

The instrumental signature present in the raw S(3839) indices, whose feature bandpass is 
3846 - 3883 A (shifted by 0.68 A), must be removed prior to a comparison with those predicted 
by models. This was carried out following Briley & Cohen (2001) by fitting the continua of the 
observed spectra with cubic splines within IRAF’s “splot” facility over the range 3600 - 5200 A. 
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The I(CH) indices determined from these continuum fits (with no CH absorption included) were 
effectively zero (—0.008±0.007) as is expected from a two-sided index. 

For S(3839), indices were computed from these continua fits and the average (0.156±0.019) 
used as a zero-point offset (as follows from the logarithmic nature of S(3839)) in our comparisons. 
Stars with spectral coverage ending redward of 3820 A were not included. The process used to 
determine the zero point offset for the uvCN indices thus assumes that the transmission T{X) of 
the atmosphere, telescope and instrument (including slit losses) normalized to the transmission at 
some fiducial wavelength, e.g. 4000 A, is a function of wavelength which does not change among 
the spectra of the M5 stars over the relevant wavelength range, i.e. that of the two bands used to 
define S(3839), from 3840 to 4000 A; this is a less stringent requirement than would be needed to 
flux the spectra. 

Slit alignment errors which vary from star to star (possible with the narrow 0.7 arcsec wide slits 
used here) combined with atmospheric dispersion are the most likely way to introduce a variation 
in the normalized T(A) from star to star, which would in turn introduce a dispersion in S(3839) 
through variation of the zero point from that adopted above. The two slitmasks used for our M71 
sample were observed at mean airmasses of 1.08 and 1.25 respectively. This should be adequate to 
avoid serious problems with atmospheric dispersion. Cohen & Cromer (1988) discuss this issue in 
more detail. 

We have tested the correlation between fit uv continuum slope and various parameters, in¬ 
cluding RA, Dec, luminosity and S(3839). The correlations range from not statistically significant 
to marginally significant, with the correlation between uv continuum slope and S(3839) being the 
most significant (correlation coefficient of 0.48) and that with luminosity being the least significant 
(correlation coefficient 0.09). 

The small corrections derived from these correlations have a maximum value of about 0.015 
and a minimum value of 0.000, which is less than 5% of the total range and smaller than the quoted 
uncertainties. Given that the S(3839) variations among the stars in our sample of subgiants in M5 
cover a range of roughly 0.3 in S(3839), these small corrections, if implemented, would have little 
affect on our resulting C and N abundances, and we have elected not to include them. 

As can be seen in Figure Q, a substantial star-to-star variation is seen among I(CH) indices 
for the stars in our M5 sample, even if one considers only subgiant branch (SGB) stars in similar 
evolutionary states. To aid in interpretation, the stars in Figure Q have been arbitrarily divided 
into CH-weak (open symbols) and strong groups (filled symbols). 

Six of the stars with the most extreme CH or CN indices are denoted as anomalous. Of these, 
five have radial velocities, all consistent with membership in M5. Four of these stars are in the 
sample at the base of the RGB and show very strong I(CH) with very weak CN. C18225_0537 has 
V = 18.06 and shows an enormous I(CH) with a strong S(3839) as well. C18206_0533 {V = 18.42) 
shows modest I(CH) but very strong uv CN. For the brighter of these stars, the anomaly noted in 
the final summed spectra is clearly present on each individual spectrum of the object, and hence it 
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is highly unlikely that these anomalies are due to CCD defects or other such problems. We plan to 
check the membership of each of these stars with higher dispersion spectra as soon as possible, but 
at the present time, we believe that most, if not all, of these anomalous stars are in fact members 
of M5. 

We proceed by restricting our attention to the subgiants at the base of the RGB in M5, 
specifically to those with 16.5 < V < 17.5 mag. Figure displays I(CH) versus S(3839) for this 
group of 43 stars. Again the very large range in each index within this small range in evolutionary 
state, hence small range in Fe// log( 5 '), is apparent. Also apparent is a strong anti-correlation 
between the strength of the CH band and that of the CN feature. The correlation is not perfect, and 
tends to turn over among the stars with the weakest I(CH). The bimodal distribution of CN band 
strengths found among the most luminous M5 giants by Smith &: Norris (1983) is also not apparent 
in the present sample. This may be the result of the higher temperatures (and correspondingly 
weaker CN bands) of our stars, combined with a finite observational error tending to blur out the 
bimodality. 

The large range in C abundances which we suspect to be present in the M5 subgiant sample 
creates an unusual situation with regard to the expected strength of the CN features. Normally, 
since there is more carbon than nitrogen, the N abundance controls the amount of CN. However, 
if C is highly depleted, there can be fewer carbon atoms per unit volume than nitrogen atoms, 
and C will control the formation of CN, as suggested by Langer (1985). Within the small range of 
Tgff covered by the M5 subgiant sample, we might expect the observed relationship between the 
molecular band indices of CH and of CN to be non-monotonic, as is in fact shown in Figure |^. In 
particular, in this figure, the stars with the weakest CH do not have the strongest CN indices. 


4.1. The Main Sequence Turn Off Stars 

There are twelve stars in our sample with V > 18.2 which are believed to be members of M5 at 
or slightly fainter than the main sequence turnoff (MSTO). Their molecular bands are, as expected, 
much weaker due to their higher Fe//- However, it does appear that variations in I(CH) are still 
present, although not as large in extent as are seen among the subgiants. Figure ^ presents the 
region of the G band of CH in four of the main sequence stars in M5 in order of increasing V mag. 
They span a total range of 1.4 mag in V from the top of the main sequence downward in luminosity 
and cooler in ^e//- The spectra have been slightly smoothed using a Gaussian with FWHM of 5 
pixels, then a continuum was fit from 3900 to 4500 A, and the spectra were normalized. The signal- 
to-noise ratio of the spectra decreases as the magnitude increases. The GH band is marginally 
stronger in the first star than in the second, and strongest in the coolest star, much stronger 
than in a main sequence star only 0.2 mag brighter. The strongest atomic features, particularly 
the Ca I lines at 4226 A and the Fe I lines at 4046 and 4384 A, are also becoming stronger as ^e// 
decreases. Higher precision data and a larger sample of main sequence turnoff stars will be required 
to determine in a definitive way the trends among CH and CN. 
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5. Comparisons with Synthetic Spectra 

Clearly the pattern of abundances underlying the CH and CN band indices of Figure ^ cannot 
be interpreted on the basis of band strengths alone - we must turn to models. The technique 
employed is similar to that of Briley & Cohen (2001), where the region of the CMD of interest is 
fit by a series of models whose parameters are taken from the Bergbusch &: VandenBerg (1992) 
0-enhanced isochrones. For M5, the 16 Gyr isochrone with [Fe/H] = -1.26 was used to select a set 
of representative model points which are listed in Table Model stellar atmospheres were then 
generated using the Marcs model atmosphere program (Gustafsson et al. 1975) at the T^ff, log(g) 
of these points. From each model, synthetic spectra were calculated using the SSG program (Bell 
& Gustafsson 1978; Gustafsson Sz Bell 1979; Bell & Gustafsson 1989; Bell, Paltoglou, & Tripicco 
1994) and the line list of Tripicco &: Bell (1995). Each synthetic spectrum was initially computed 
from 3,000 to 12,000 A in 0.1 A intervals and convolved with V and I filter curves as described in 
Gustafsson &: Bell (1979) and Bell &: Gustafsson (1989) to yield colors appropriate to each model. 

Figure ^ shows these model points superposed on the M5 V, I CMD. We assume stellar masses 
as given by the isochrone, a microturbulent velocity of 2 km/s, and the cluster parameters given in 
§|I|. As is apparent from the Figure, the isochrone closely follows the M5 CMD. We have compared 
Teff derived from the isochrone itself, which values we subsequently adopt, with T^fj derived from 
the observed V — I colors of each of the subgiants individually. We use the grid of predicted 
colors from Houdashelt, Bell & Sweigart (2000) for this purpose, together with the log(( 7 ) from the 
isochrone. (The dependence on log( 5 f) is relatively small.) The mean difference is ~ +70K, with a 
large a of 90 K. While some part of this mean difference might be photometric errors, if interpreted 
as an error in matching the color of the isochrone at a hxed V with the locus of the M5 stars, this 
corresponds to a systematic difference inV — I of about 0.03 mag. A hint of this is seen in Figure]^, 
where the cluster locus does appear slightly redder than the adopted model points near F~17.5. 
Ignoring these small differences, we therefore believe that assigning Tgff,log{g) utilizing the V mag 
plus the isochrone is a reasonable way to represent the mean properties of the M5 subgiant sample 
as a function of luminosity. 

Using these models we have further calculated a series of synthetic spectra at higher resolution 
(0.05 A intervals) from 3500 to 5500 A with a variety of G, N, and O abundances. These spectra were 
then smoothed to the resolution of the observed spectra and the corresponding I(GH) and S(3839) 
indices measured. Zero point shifts were also measured from the synthetic spectra by loading 
them into IRAF and performing the continuum fits discussed above for the observed spectra. The 
resulting offset for S(3839) was 0.046±0.010 and 0.003±0.001 for I(GH). 

The values of I(GH) from five sets of models with [G/Fe] from -0.35 to -1.10, [0/Fe] = +0.25, 
and G^^/G^^= 10 are plotted with the observed indices in Figure |^. The spread in [G/Fe] among 
the M5 SGB stars appears well represented by this range, which is very similar to that reported 
among a sample of more luminous M5 SGB stars by Briley et al. (1992) (although a lower value 
of [Fe/H] = -1.40 was chosen in their study). We note the very low G abundances implied by the 
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weakest I(CH) indices among the SGB sample, which are consistent with the low C abundances 
observed among the brightest M5 giants (Smith et al. 1997). However, in the albeit small Smith 
et al. sample of five giants, no analog to the CH-strong stars (at least in terms of [C/Fe]) were 
observed. 

The large range in C abundances indicated in Figure ^ leads to a difficulty in assessing N 
abundances via CN band strengths, as discussed in §^. Among the CH-strong SGB stars (see the 
left panel of Figure P) there is a small (~ 0.1 dex) spread in [C/Fe], with N being the minority 
species controlling the formation of CN. Thus, the CN band strengths of these stars can be assumed 
to roughly represent the [N/Fe] abundances. This is demonstrated in Figure ^ (left), where the 
S(3839) indices are compared to those from two series of synthetic spectra, the first with [C/Fe] = 
-0.35, [N/Fe] = -1-0.10 and the second with [C/Fe] = -0.45, [N/Fe] = -1-1.00 (with [0/Fe] = -1-0.25 
and C^^/C^^=10 for both, as in Figure ^). Clearly the star-to-star spread in N abundances of the 
CH-strong SGB stars approaches a factor of 10. 

For the CH-weak stars, the situation is more complex - the lower C abundances can regulate 
CN formation. Thus, as suggested by the discussion of the large spread in [C/Fe] among the CH- 
weak stars essentially destroys any relation between S(3839) and N abundance. This is illustrated 
in Figure ^ (right) where S(3839) indices from model spectra with three different combinations of 
[C/Fe] and [N/Fe] are plotted. While [C/Fe] = -0.65, [N/Fe] = -1-1.20 results in the stronger of the 
CN bands, the models with [C/Fe] = -0.55, [N/Fe] = -1-0.75 have CN band strengths essentially 
indistinguishable from [C/Fe] = -1.10, [N/Fe] = -1-1.30. Indeed, for the most C-poor of the current 
SGB stars, significant enhancements of N are required to reproduce the observed S(3839) indices. 

We defer an attempt to deduce C and N abundances among the M5 main sequence turn off 
stars until a larger sample of spectra of such stars with higher signal-to-noise ratio than the present 
set of spectra becomes available. 


6. Inferred C and N Abundances Among the Subgiants 

To disentangle the underlying C and N abundances from the CH and CN band strengths, we 
have simultaneously fit the [C/Fe] and [N/Fe] abundances corresponding to the observed I(CH) 
and S(3839) indices of the SGB stars. The temperatures and gravities used for the models were 
chosen based on the V magnitudes of the observed stars in the present photometry as fit to the 
isochrone in Figure ^ using cubic splines. The validity of this procedure was established in §^. The 
abundances of C and N were then varied until the model I(CH) and S(3839) indices matched the 
observed values as closely as possible, including the derived zero points. The mean rms error in 
the fits to the observed indices was 0.005. For all stars, a value of [0/Fe] = -1-0.25, C^^/C^^= 10, 
and a microturbulent velocity of 2 km/s was assumed (the sensitivity to these assumptions will be 
explored below). The resulting C and N abundances for the M5 SGB stars are listed in Table ^ 
and plotted in Figure 0 
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Immediately apparent in Figure |T^ is a dramatic anti-correlation between the C and N abun¬ 
dances of the M5 SGB stars. Also plotted are the [C/Fe] and [N/Fe] abundances of a sample of 
more luminous M5 SGB stars reported in Briley et al. (1992) (note we have plotted the abundances 
from their [Fe/H] = -1.25 models), which compare very well with the present results. The error 
bars were determined by repeating the fitting process while including shifts in the observed indices 
of twice the average error among the SGB indices as derived from Poisson statistics (0.005 in I(CH) 
and 0.02 in S(3839)). The shifts were included in opposing directions (e.g., -1-0.005 in I(GH) and 
-0.02 in S(3839), followed by -0.005 in I(GH) and -1-0.02 in S(3839)) and reflect likely errors in the 
abundances due to noise in the spectra. 

Of much greater concern are the systematic errors, as they have the potential to exaggerate 
any C versus N anticorrelation (i.e., stars with overly low C abundances will naturally require 
greater abundances of N to match the observed GN band strengths). To assess the role of many 
of our assumptions in shaping Figure 0, we have chosen three representative SGB stars and again 
repeated the fitting of the CN and GH band strengths with different values of [Fe/H], [0/Fe], 
G12 /c 13^ etc. These results are presented in Table where it may be seen that the sensitivity of 
the derived G and N abundances to the choice of model parameters is remarkably small (well under 
0.2 dex for reasonably chosen values), as would be expected from these weak molecular features. 
We have also plotted the SGB G and N abundances as both functions of V and V — I colors in 
Figure |l^ to evaluate possible systematic effects with luminosity (albeit over a small range) and 
temperature; none appear to be present. 

Evaluating the absolute abundance scale is more difficult as external comparisons are limited. 
For the main sequence stars in 47 Tuc, we can compare the results of Briley et al. (1991, 1994), car¬ 
ried out in a manner fairly similar to the present work, with the independent analysis of a different 
sample of stars by Gannon et al. (1998). This suggests we may be systematically underestimating 
the absolute G abundance by about 0.15 dex, and overestimating the N abundance by about 0.2 
dex. However, such a shift cannot account for the large range in C and N abundances apparent 
in Figure H- We therefore conclude that the G versus N anti-correlation among the SGB stars in 
Figure 0 is indeed real. 


7. Review of Previous Studies of C and N in Globular Clusters 

7.1. The Case of M5 

There have been several previous studies of the C and N abundances in M5. Smith & Norris 
(1983) measured CH and GN indices for a sample of 29 stars near the RGB tip. They found a 
bimodal distribution of GN absorption. Langer et al. (1985) reobserved six stars selected from this 
sample, three CN-strong and three GN-weak, and calculated the change in mean C and N abundance 
between the two groups. Briley et al. (1992) observed 14 subgiants with 16.2 <V< 16.5, analyzed 
in a manner similar to the present paper, to find large stochastic variations from star-to-star of 
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a magnitude similar to those found here for still fainter subgiants, and with a clear detection of 
bimodality for the S(3839) indices of these subgiants. 

A question of considerable import for understanding the behavior of C and N in M5 is whether 
or not there is any change in the mean C abundance as one moves from the RGB tip down to the 
base of the subgiant branch. Any effect dependent on evolutionary state is an argument favoring 
internal mixing rather than an external (primordial or pollution) origin for this phenomenon. While 
this issue is complicated by the need to convert between three different definitions for the strength 
of absorption at the G band of CH, the best that can be said at the present time is that any 
difference in the mean [C/H] between the tip and the base of the RGB in M5 must be less than 0.3 
dex. We seem to be in the peculiar position of having a better sample of stars with published C 
and N abundances between 16.0 < V < 17.5 than exists at the RGB tip. We have already initiated 
an effort to remedy this and to refine this crucial number. 


7.2. Other Globular Clusters 

Gohen (1999) has shown that the C/N anti-correlation seen among bright RGB stars in M71 
extends to the stars at the main sequence turnoff and even fainter in M71. At the level of the main 
sequence, the distribution of both the CH and CN indices appears to be bimodal. Briley & Cohen 
(2001) have used model atmospheres combined with synthetic spectra in an approach similar to 
that of the present paper to show that the C and N abundance range seen at the level of the main 
sequence is comparable to that seen among the bright red giants of M71 by many previous studies, 
the most recent of which is Briley, Smith & Claver (2001). To reproduce the M71 CH and CN 
bands requires a range in C of a factor of 2 starting from [C/Fe] = 0.0 and a range in N of a factor 
of 10, starting from an initial N enhancement of a factor of 2.5 ([N/Fe] = -1-0.4 dex, together with 
a decrease in O from [0/Fe] = -1-0.4 to -1-0.1. These numbers are very similar to what is required 
to explain the M5 subgiants, with the exception of [C/Fe], which exhibits a larger variation among 
the M5 subgiants (a factor of 6 versus 2). 

Norris, Freeman & DaCosta (1984) studied the CH and CN indices among a sample of 112 
bright RGB stars in 47 Tuc (a high metallicity cluster, with [Fe/H] similar to that of M71), while 
Cannon et al. (1998) observed a sample of comparable size from the base of the RGB to the upper 
main sequence. Large stochastic star-to-star variations with similar anti-correlations between C 
and N abundances and with a bimodal behavior of CN were found at all luminosities examined. 
The range of mean abundances between the CN-strong and CN-weak groups appeared unchanged 
from the lower RGB to the main sequence, with the CN-strong group having <[C/H]>= —0.15, 
<[N/H]>= -1-1.05 dex, as compared to <[C/H]> = -1-0.06, <[N/H]>= -1-0.20 dex for the CN-weak 
stars. The N range is comparable to that found here for M5, but again, as with M71, the C range 
is a little smaller in 47 Tuc. Similar to the situation in M5, large samples exist both at the tip of 
the RGB and at the main sequence in 47 Tuc, but they are not well tied together, and it is not 
possible from the published papers themselves to establish the magnitude of any systematic change 
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in the mean < [C/H] > with luminosity, beyond a statement that it cannot be large. 

Determinations of the ratio have been made for the brightest RGB stars in 47 Tuc 

by Brown, Wallerstein & Oke (1990) and Bell, Briley, & Smith (1990), as well as in M71 by Briley 
et al. (1997). They obtain ~4 - 8, correlated with CN-band strength, implying that the 

surface C-poor, N-rich envelope material has been exposed to proton capture. 

The seminal survey of C and N in the intermediate metallicity clusters M3 and M13 (two 
clusters of slightly lower metallicity than M5) by Suntzeff (1981) established the existence of a 
bimodal distribution of CN band strength, at least for My ^ —0.4, with stochastic variations in 
C and N abundance from star to star in each cluster. Smith et al. (1996) combine their C and 
N abundances for a sample of stars near the tip of the RGB in each cluster with published O 
abundances to show that the total of G+N+0 is constant for these luminous giants. 

The only other cluster with similar data available for its MSTO stars is NGG 6752, a cluster of 
similar metallicity to M3 and M13. Suntzeff & Smith (1991) observed both GH and CN variations 
among its MSTO stars. Gratton et al. (2001) report an anticorrelation between O and Na, as well 
as between Mg, Al and C and N, similar to the relations found among the luminous giants. Suntzeff 
& Smith (1991) also found low (3 - 10) ratios, which they attributed to some mixing taking 

place but with primordial variations already in place. 

M92 was studied in detail by Carbon et al. (1982), who used the NH band at 3360 A, as 
the CN bands become too weak at such low metallicities. They studied 45 giants and subgiants 
with My < +2 and found strong stochastic star-to-star variations over their full luminosity range of 
about a factor of 3 in C and a factor of 10 in N. However, C and N are not in general anti-correlated. 
The data for fainter stars is limited at present (at least by the standards of the large samples we 
have achieved for M71 and for M5). However, it is clear that in addition to the variations, there 
is an easily detectable systematic decrease in the mean C abundance by about a factor of 10 from 
the base of the RGB to the top of the RGB (hanger et al. 1986, Bellman et al. 2001). They argue 
that this is the key fact in establishing that internal mixing is the dominant effect controlling the 
C and N abundances in M92. 

M15, with metallicity similar to M92, was studied by Trefzger et al. (1983). In a sample 
of 33 bright giants reaching to My = +1.2, they found strong stochastic star-to-star variations, 
and found that the mean <C/Co > declines with advancing evolutionary state, but the mean N 
abundance does not change. Furthermore, they found that N/Nq was so large for about 1/3 of the 
sample that it could not be explained even by converting all C into N. 

GH in the upper RGB stars in NGG 6397, another metal poor cluster, was explored by Bell, 
Dickens &: Gustafsson (1979). 
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8. Discussion 


The primary facts that we have established for M5 are that there are strong stochastic varia¬ 
tions from star-to-star of both C and N, with C and N anti-correlated, similar to those seen in M71. 
These variations are definitely present among the subgiants at the base of the RGB to My^ + 3 
and appear to extend to the main sequence stars as well. A question of considerable import for 
understanding the behavior of C and N in M5 is whether or not there is any change in the mean 
C abundance as one moves from the RGB tip down to the base of the subgiant branch such as 


is clearly seen in M15 and M92. As discussed in §7.2 above, the best that can be said at the 
present time is that, unlike the case of M92 or M15, any systematic decline in the G abundance 
with evolutionary state as one moves up the RGB appears to be small, < 0.3 dex. 


8.1. Implications for Stellar Evolution 

A classical review of post-main sequence stellar evolution can be found in Iben & Renzini 
(1983). Their description of the consequences of the first dredge up phase, the only dredge up 
phase to occur prior to the He flash, indicates that a doubling of the surface and a 30% 
reduction in the surface can be expected, together with a drop in the ratio of from 

the solar value of 89 to ~20, as well as a drop in surface Li and B by several orders of magnitude. 
Observations of field stars over a wide range of luminosities conform fairly well to this picture, see 
e.g. Shetrone et al. (1993), Gratton et al. (2001), although additional mixing of Li and lower than 
predicted ratios of seem to occur even among field stars (do Nascimento et al. 2000). 

To match the observations of variations in abundances among globular cluster red giants which 
far exceed those described above, additional physics must be introduced into calculations of dredge 
up in old metal poor stars. Relevant phenomena include meridional mixing as described by Sweigart 
&: Mengel (1979) as well as turbulent diffusion (see Charbonnel 1994, 1995) and the insights of 
Denissenkov &: Denissenkova (1990) concerning the importance of the ^^Ne(p, 7 )^^Na reaction as a 
way to produce p-burning nuclei. 

The clear prediction of the most current calculations of this type by Denissenkov &: Weiss 
(1996), Cavallo, Sweigart & Bell (1998) and Weiss, Denissenkov & Gharbonnel (2000) is that the 
earliest that deep mixing can begin is at the location of the bump in the luminosity function of 
the RGB which occurs when the H-burning shell crosses a sharp molecular weight discontinuity. 
Zoccali et al. (1999) have shown that the luminosity of the RGB bump as a function of metallicity 
as determined from observation agrees well with that predicted by the theory of stellar evolution. 
Bono et al. (2001) further suggest that the agreement between the predicted luminosity function 
and actual star counts along the RGB in the vicinity of the bump in a suite of globular clusters 
is so good that mixing cannot have occurred any earlier, otherwise the evolutionary lifetimes, and 
hence the observed luminosity function, of such stars would have been affected by the mixing of 
He. 
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Zoccali et al. (1999) give the expected location of the RGB bump in M5 to be 0.3 mag brighter 
than the HB, i.e. at R ~ 14.8 or My ~ +0.5. Yet we see strong star-to-star variations in C and 
N abundances as well as a strong anti-correlation between them within a large group of cluster 
members at R ~ 17.1 (My ~ +2.8), more than 2.3 mags fainter at the base of the RGB. We see 
hints of such variation continuing on the upper main sequence at My ~ +3.7. 

The range of luminosity over which these C and N variations are seen is becoming more and 
more of a problem for any scenario which invokes dredge up and mixing. Unless we have missed 
some important aspect of stellar evolution with impact on mixing and dredge up, we must declare 
the mixing scenario a failure for the specific case of M5 from our present work and M71 from 
our previous work (and several other globular clusters as well from the work of others). Even 
the theoreticians in the forefront of this field are beginning to admit that deep mixing alone is 
not sufficient (Denissenkov &: Weiss 2001, Ventura et al. 2001). Unless and until some major new 
concept relevant to this issue appears, we must now regard the fundamental origin of the star-to-star 
variations we see in M5 as arising outside the stars whose spectra we have studied here. 

The strong anti-correlation between C and N, however, does suggest that GN-cycle material 
must be involved, and that this material has somehow reached the surface of these subgiant stars 
in M5. Since we know it cannot come from inside these stars, it must come from some external 
source. As reviewed by Lattanzio, Gharbonnel &: Forestini (1999), CN and ON cycling is known 
to occur in AGB stars, and AGB stars are also known to have sufficient dredge up to bring such 
material to their surfaces. We might speculate that the site of the proton exposure could be a 
previous generation of high mass stars, which then suffered extensive mass loss (either in or outside 
of binary systems) and polluted the generation of lower-mass stars we currently observe, while the 
higher mass stars are now defunct. 


8.2. ON Burning 


Let us adopt as a working hypothesis that the C and N abundance variations we are seeing in 
the present subgiants are the result of the incorporation of material exposed to the GN-cycle (i.e., 
proton capture reactions) in a now evolved population of more massive (2 - 5Mq) AGB stars, as was 
originally suggested in D’Antona, Chieffi &: Gratton (1983). Indeed, recent models of metal-poor 
AGB stars by Ventura et al. (2001) suggest temperatures at the bases of the convective envelopes 
of such stars are capable of these reactions. The new generation of precision abundance analyses 
of globular cluster stars over a wide range of luminosity such as that of Ramfrez Gohen (2002) 
for M71 demonstrate that the abundances of the a-capture and s-process elements are constant, 
and so we further assume that the early cluster environment of M5 was not significantly polluted 
by the ejecta of even more massive stars. 


We now investigate whether this hypothesis is consistent with the C and N abundances we have 
derived for the main subgiant sample of M5. Figure 12 shows the sum of the C and N abundance 
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as a function of the C abundance of the sample of M5 subgiants. The solid dot shows the predicted 
location assuming the initial C and N abundances (Cq, Nq) are the Solar values reduced by the 
metallicity of M5 ([Fe/H] = —1.2 dex). Thus this is the initial location for no burning and for 
a Solar C/N ratio. If the present stars incorporated material in which just C was burned into 
N, then the locus of the observed points representing the M5 subgiant sample should consist of a 
single horizontal line, with the initial point, the presence of no CN-cycle exposed material, at the 
right end of the line (the maximum C abundance) and the left end of the line corresponding to a 
substantial fraction of the star’s mass (i.e. the atmosphere plus surface convection zone) including 
C-poor, N-rich AGB stellar ejecta. Furthermore, if the initial C/N ratio of the cluster is not Solar, 
then the locus should still be a horizontal line, but located at a different vertical height in this 
figure. 

The maximum possible N enhancement for a cluster SGB star with these assumptions occurs if 
the star formed entirely from AGB ejecta in which all C has been converted into N. For initial values 
(Co, Nq) (not expressed as logarithms), this maximum N enhancement would be (Cq + No)/No. 
If the initial value was the Solar ratio, Cq/Nq ~ 3.2, the resulting maximum N enhancement is a 
factor of ~4.2, while for an initial Cq/Nq of 10, the maximum N enhancement is a factor of 11. 

Now we examine the behavior of the C and N abundances among the M5 subgiant sample as 
inferred from our observations. It is clear that the assumption that the only thing happening is 
inclusion of material in which C was burned into N must be incorrect. The sum of C+N seems 
to systematically increase by a factor of ~5 between the most C rich star and most C deficient 
star. The discussion of the errors, both internal and systematic, in suggest maximum systematic 
errors of —0.2 dex for log(C/H) and +0.2 for log(N/H). This is completely insufficient to explain 
such a large trend as errors. 

We thus have a serious discrepancy. The sum of C+N was not constant as C was burned in 
the AGB sites. Furthermore the observed range in N abundances is very large. The most obvious 
way to reproduce this is to include O burning as well as C burning. If we adopt Solar ratios as our 
initial values, then a substantial amount of O burning is required. 

Figure suggests that the initial ratio of C/N is not quite Solar, although not too far off. 
Adopting the Solar value as the initial C/N ratio, we calculate the minimum amount of O which 
must be burned at the base of the AGB envelopes to reproduce the locus observed in the figure 
(under the arguable assumption of the most extreme of our stars having formed largely from such 
material - this will, however, provide us with at least an estimate of the minimum burning required). 
We need to produce a N enhancement of a factor of 10. The Solar ratio is C/N/0 = 3.2/1/7.6, 
so if all the C and 50% of the O were converted, we have an enhancement of N of a factor of 8 
available to the present stars. Oxygen is typically found to be overabundant with respect to Fe 
in old metal-poor systems (see Melendez, Barbuy & Spite 2001, Gratton et al. 2001, Ramirez & 
Cohen 2002, and references therein); we assume [0/Fe] ~ +0.3 dex, a typical value. Then the 
initial C/N/0 ratios will be 3.2/1/15.2. Note that the same amount of O has to be burned to 
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produce the observed distribution of C and N abundances, but in this case it is a considerably 
smaller fraction of the initial O. 

Returning to the AGB models of Ventura et al. (2001), the requirement for substantial O 
burning that emerges from our analysis of the CH and CN bands in the M5 subgiants may not 
be an unreasonable constraint - for metal-poor AGB stars they find temperatures sufficient for 
CNO-processing at the bases of AGB stars in a wide mass range. Following their Z=0.001 ([Fe/H] 
= -1.3) models we find surface O abundances dropping by a factor of 2 to 20 in masses from 4.5 
to 2Mq. We also note that under the assumption of little change in [C/Fe] (less than 0.3 dex 
as discussed above) taking place during the RGB ascent of the present low mass M5 stars, one 
should also expect little change in [0/Fe] as well, and that the O abundances of the present bright 
giants reflect their “primordial” values. The observed [0/Fe] abundances of the bright M5 giants 
by Sneden et al. (1992) are not inconsistent with this idea. Their “0-rich” stars average around 
[0/Fe] = -1-0.3 while their most “0-poor” stars are depleted by a factor of 3.5. If this is the result 
of ON-cycle exposure, more than enough N can be produced to explain the present results. 

This simple test is of course leaves significant questions unanswered. Problems include whether 
there were a sufficient number of AGB stars present to return the required quantity of material 
and the efficiency of any mechanism to incorporate it in the present stars. This is a non-trivial 
issue if the most C/O-poor, N-rich SGB stars formed with a preponderance of AGB ejecta - to 
reduce the C abundance by a factor of 6 by adding C-poor ejecta would require some 83% of 
the present star’s mass to be made of this material. Note, as has been pointed out by several 
authors, these abundance inhomogeneities cannot simply be surface contaminations as they would 
be diluted by the increasing depth of the convective envelope during RGB ascent. Also, the range 
of C abundances among the M5 SGB stars appears much larger than that of any cluster studied 
to date. While this can perhaps be explained with regard to the more metal-rich clusters, whose 
polluting AGB stars should have undergone less ON-cycling, other even more metal-poor clusters 
than M5 appear to have smaller star-to-star C variations among their less evolved stars. This can 
be seen in the [C/Fe] versus luminosity diagrams of Carbon et al. (1982) and Bellman et al. (2001) 
for M92, where the range in [C/Fe] among the least luminous stars is relatively small (certainly not 
the factor of 6 seen here in M5). But this may be perhaps explained by varying the efficiency for 
incorporating AGB ejecta into subsequent generations of stars among the proto-globular clusters. 

Following Ventura et al. (2001), we predict substantial O variations, anti-correlated with Na 
and Al, to be present among the M5 subgiants and fainter stars. If the source of the proton 
exposed materials is indeed moderate mass AGB stars, these inhomogeneities should also correlate 
with Li abundances among the less evolved stars. We expect the verification or lack thereof of these 
predictions to be available shortly. 

Another important point is that Figure ^ shows no evidence for any bimodality in the distri¬ 
bution of the C and N abundances for the small section of the subgiant branch in M5 covered by our 
sample. The distribution along the locus of abundance appears to first order to be uniform. There 
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is no preponderance for stars populating the extremes of high C or low C. An artificial suggestion 
of bimodality, or more correctly a tendency towards high CN strengths, could be produced by the 
distribution of C and N abundances shown in Figure 0 or by saturation effects in the molecular 
bands themselves, which would only become apparent in cooler stars with stronger molecular bands. 

We also recall the anomalous stars in our sample, most of which are believed to be members 
of M5. Two of these in particular have enormously strong CH bands; we can offer no explanation 
for these stars at present. 


8.3. Additional Implications for Stellar Evolution 

Figure ^ was used above to demonstrate that ON burning is required by considering the 
required N enhancement factor. This figure also displays a correlation between the sum of C and 
N number densities with the C abundance, i.e. the C abundance is correlated with the C/N ratio 
among the M5 subgiant stars. Here we explore the consequences of the existence of this correlation 
for the origin of the C and N variations themselves. Any external mechanism for producing these 
variations will involve an efficiency factor for the incorporation of material. We expect this factor 
to depend on the mass of the star itself, how much additional mass is incorporated (AM), and 
the initial C and N abundances in the star itself and within AM. Since these properties of AM 
might be expected to fluctuate wildly, this process therefore should show a lot of stochastic random 
variability. 

It is easy to imagine that various parcels of AM have a wide range in C/N ratios due to 
the varying amount of nuclear processing that each might have experienced, and thus the strong 
stochastic star-to-star variations in C and N as well as their overall anti-correlation that we observe 
among the M5 subgiants can be reproduced. However, the correlation above requires a correlation 
between the mass of the parcel accreted, AM, and the C/N ratio of the material within this mass, 
and that seems to be rather artificial given the random nature of the process. We thus conclude that 
is difficult to reproduce the correlation described above with such an external mechanism involving 
accretion of “polluted” material from AGB stars. Unless we have made gross errors in the C and 
N abundances in our M5 subgiant sample far beyond what we believe might have occured, the 
existence of this correlation suggests that we should not totally rule out internal mechanisms as 
yet. 


9. Summary 

We have presented photometry and spectroscopy for a large sample of M5 subgiants and several 
main sequence turn off stars. An analysis of the SGB spectra reveals significant and anti-correlated 
star-to-star variations in C and N abundances, as would be expected from the presence of proton- 
capture exposed material in our sample stars. Similar variations in CH also appear to be present 


- 21 - 


in the main sequence turn off spectra, but the signal in the current sample is too low for a detailed 
analysis. The evolutionary states of these stars are such that the currently proposed mechanisms 
for in situ modifications of C, N, O, etc. have yet to take place. On this basis, we infer that the 
source of proton exposure lies not within the present stars, but more likely in a population of more 
massive (2 - 5Mq) stars which have “polluted” our sample. 

The C+N abundances derived for the M5 subgiants show a large systematic increase in C+N 
as C decreases. To reproduce this requires the incorporation not only of CN, but of ON-processed 
material as well. Furthermore, the existence of this correlation is quite difficult to reproduce with 
an external mechanism such as “pollution” with material processed in a more massive AGB star, 
which mechanism is fundamentally stochastic in nature. We therefore suggest that although the 
internal mixing hypothesis has serious flaws, new theoretical insights are needed and it should not 
be ruled out yet. 


The entire Keck/HIRES and LRIS user communities owes a huge debt to Jerry Nelson, Gerry 
Smith, Steve Vogt, Bev Oke, and many other people who have worked to make the Keck Tele¬ 
scope and HIRES and LRIS a reality and to operate and maintain the Keck Observatory. We are 
grateful to the W. M. Keck Foundation for the vision to fund the construction of the W. M. Keck 
Observatory. The authors wish to extend special thanks to those of Hawaiian ancestry on whose 
sacred mountain we are privileged to be guests. Without their generous hospitality, none of the 
observations presented herein would have been possible. 
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Fig. 1.— The histogram of the measured heliocentric radial velocities for approximately half of the 
sample stars in M5 is shown. The cluster Vr of Pryor & Meylan (1993) (+51.9 km s“^) is indicated 
as a dashed vertical line. 
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Fig. 3.— Part of the LRIS-B spectra of two probable members of M5. The stars are essentially 
identical in V mag and V — I colors and are both located at the base of the subgiant branch. Note 
the difference in the strength of the G band of CH in these two spectra. 
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Fig. 4.— The measured I(CH) and S(3839) indices are plotted for the program stars as a function 
of V. The sample has been arbitrarily divided into two gronps: CH-strong (filled markers) and 
CH-weak (open). Large and significant star-to-star differences exist in both CH and CN band 
strengths among the SGB stars. The decreasing spread in indices with luminosity is the result of 
increasing temperatures near the main sequence turn off. The error bar shown in the upper-right 
box represents the uncertainty in determining the slope of the uv continuum as described in §^. 
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Fig. 5.— The CH and uvCN band indices are plotted against each other for the 43 stars with 
16.5 < V < 17.5 mag at the base of the RGB in M5. The two stars with the strongest CH indices 
were classified as anomalous. 
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Fig. 6.— The region of the G band of CH in four main sequence stars in M5 going from top to 
bottom in order of increasing V mag. A constant arbitrary vertical offset between the spectra has 
been imposed for clarity. The label for each object gives in parentheses the DN/pixel at the blue 
and red end respectively in the original spectrum so that the signal-to-noise ratio can be calculated. 
The spectra have been slightly smoothed, then a continuum was fit from 3900 to 4500 A, and the 
spectra were normalized. In addition to looking at the strength of the CH band, note that the 
strongest atomic features, particularly the Ca I lines at 4226A and the Fe I lines at 4046 and 
4384 A, are also becoming stronger as Tg// decreases. 





- 32 - 


CD 



0 0.5 1 1.5 

V—I (mag) 


Fig. 7.— The model points for the grid of stellar atmospheres and spectral syntheses obtained 
from the isochrone are shown superposed on the (F, I) CMD diagram of M5 in the region of the 
subgiant branch and the upper main sequence. 
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Fig. 8.— The observed I(CH) indices are plotted with calculated band strengths for several differing 
C (and N) abundances as listed in Table ^ as a function of V magnitude. As in Figure the stars 
have been divided into two groups based on CH-band strength. Note that the range of CN-band 
strengths observed requires nearly a 0.75 dex star-to-star variation in [C/Fe] among the SGB stars. 
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Fig. 9.— Observed S(3839) indices are plotted with the model indices from Table Q against V 
magnitude. As with Figure ^ a significant spread in abundances is present to within 0.5 mags of 
the MSTO. Among the MSTO stars themselves, the higher temperatures have resulted in CN band 
strengths too weak to measure accurately in the present spectra. The error bar shown in the boxes 
indicates the uncertainty in the S(3839) offset due to the slope of the uv continuum (see §^. 
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Fig. 10.— The resulting [C/Fe] and [N/Fe] abundances for the M5 SGB stars in Table ^ are plotted. 
A strong C versus N anti-correlation is evident which also compares well with the relation from 
Briley et al. (1992) among a sample of more luminous cluster members (plotted as filled squares). 
The presence of such an anticorrelation, although suggestive of the presence of atmospheric material 
exposed to the CN-cycle, is difficult to explain via internal processes given the evolutionary state 
of the present sample of stars. 
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Fig. 11.— The derived C and N abundances are plotted against photometry for the SGB stars 
of Table |^. No systematic trends with either luminosity or temperature are apparent in the abun¬ 
dances. 
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Fig. 12.— The sum of the derived C and N abundances is plotted as a function of the C abundance. 
The large hlled circle marks the location for both C and N depleted by a factor of 16, adopting the 
abundance of M5 of [Fe/H] = —1.2 dex, with C/N at the Solar ratio. The horizontal line extending 
to the left of that represents the locus of points for C gradually being converted into N, with the 
left end of the line having C/Cq = 0.1. 




Table 1. 


Probable Non-Members of M5 


ID 

Sample 
(P or S) 

Vr 

(km s“^) 

Abs. Spectrum 

A[U,/]'^ 

(mag) 

Non-members: 

C18262_0726 

P 

-73.1 

Metal rich 

+om 

C18465_0208 

P 

b 

Metal rich 

+0.21 

Probable Members; 

C18285_0147 

P 

-h60.7 

OK 

-F0.14 

Cl8408.0805 

P 

b 

OK 

-0.14 

018488.0120 

P 

b 

OK 

-0.14 

Uncertain: 

018383.0722= 

S 

b 

BHB ? 

-0.16 


^This is \V — I for the star] ~ \y — I for the cluster locus] evaluated at 
the V mag of the star. 

'^No radial velocity is available for this star. 

'^This faint star is possibly a blue straggler in M5 or possibly a back¬ 
ground early-type star. 
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Table 2. Photometry for Members of M5 in Our Sample 


ID'^ 

B 

V 

I 


(mag) 

(mag) 

(mag 


C18149_0401 17.66 16.93 16.07 

C18159_0608 17.79 17.10 16.25 

018174.0329 15.80 14.95 13.96 

018175.0332 17.63 16.93 16.07 

018177.0109 17.76 17.10 16.20 

018177.0118 20.00 19.56 18.81 

018177.0123 18.04 17.37 16.48 

018181.0205 17.89 17.24 16.34 

018182.0207 19.17 18.81 18.11 

018184.0442 17.68 17.00 16.14 

018187.0216 17.94 17.31 16.41 

018191.0554 17.81 17.12 16.27 

018195.0301 17.77 17.07 16.20 

018200.0241 17.92 17.21 16.32 

018200.0251 18.97 18.54 17.87 

018200.0351 17.84 17.15 16.30 

018204.0418 17.52 16.83 15.97 

018204.0524 17.78 17.09 16.24 

018204.0521 18.34 17.68 16.88 

018206.0733 18.88 18.42 17.82 

018211.0559 18.59 18.06 17.42 

018225.0533 20.41 19.83 19.08 

018225.0537 17.79 17.07 16.20 

018240.0320 17.54 16.92 16.06 

018240.0621 17.65 16.93 16.09 

018240.0705 17.71 17.03 16.17 

018241.0216 17.63 17.00 16.05 

018243.0634 17.78 17.16 16.26 

018243.0638 18.87 18.38 17.73 

018246.0716 17.72 17.03 16.16 

018257.0756 18.00 17.31 16.47 

018262.0723 19.19 18.70 18.07 

018285.0147 17.82 17.32 16.24 

018268.0132 17.73 17.07 16.18 

018384.0725 17.50 16.78 15.95 
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Table 2—Continued 


ID^ 

B 

(mag) 

V 

(mag) 

/ 

(mag 

Cl8384.0728 

18.58 

18.14 

17.46 

C18386_0709 

17.60 

16.92 

16.03 

Cl8408.0805 

17.62 

16.75 

16.04 

C18412.0233 

17.78 

17.17 

16.20 

C18413.0647 

17.86 

17.26 

16.30 

C18422.0256 

18.89 

18.52 

17.82 

C18422.0306 

17.86 

17.24 

16.31 

C18422.0318 

19.73 

(19.23)^ 

18.54 

Cl8422.0748 

18.58 

18.14 

17.46 

C18422.0757 

17.74 

17.03 

16.16 

C18429.0337 

18.00 

17.40 

16.43 

Cl8437.0548 

19.13 

18.62 

17.96 

Cl8442.0658 

17.93 

17.25 

16.39 

C18448.0557 

17.75 

17.04 

16.17 

C18461.0520 

17.96 

17.32 

16.43 

C18478.0505 

18.75 

18.32 

17.67 

Cl8478.0508 

17.79 

17.10 

16.25 

Cl8484.0608 

17.82 

17.14 

16.28 

C18488.0120 

17.79 

16.85 

16.14 

C18496.0637 

17.94 

17.29 

16.42 

C18501.0410 

19.38 

18.73 

18.18 

C18502.0146 

17.65 

17.00 

16.07 

C18502.0405 

17.57 

16.91 

15.99 

Cl8502.0447 

17.68 

17.01 

16.13 

BHB Star 

C18386.0713 

15.56 

15.60 

15.59 


“^The star names are derived from their 
J2000 coordinates. Star C12345_5432 has 
coordinates 15 12 34.5 +2 54 32 (J2000). 

^The V mag is obtained from the B,B — 
I assuming the object is on the main se¬ 
quence of M5, as suggested by the available 
B,I photometry. 
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Table 3. Indices, Model Parameters, and Resulting Abundances 


Star 

Oomment 

I(GH) 

€l(GH) 

S(3839) 

eS(3839)'^ 

Te// 

log g 

[G/Fe] 

C18149_0401 

SGB-OH-Strong 

0.192 

0.004 

0.236 

0.012 

5192 

3.23 

-0.35 

C18159_0608 

SGB-OH-Weak 

0.157 

0.004 

0.229 

0.012 

5216 

3.31 

-0.52 

C18174_0329 

SGB-GH-Strong 

0.221 

0.004 

0.346 

0.012 

4849 

2.29 

-0.42 

C18175_0332 

SGB-GH-Weak 

0.131 

0.004 

0.273 

0.010 

5192 

3.23 

-0.69 

C18177_0109 

SGB-GH-Weak 

0.152 

0.004 

0.191 

0.012 

5216 

3.31 

-0.55 

C18177_0118 

MS-GH-Weak 

0.022 

0.012 

0.049 

0.042 

- 

- 

- 

C18177_0123 

Anomalous 

0.193 

0.004 

0.015 

0.016 

5260 

3.44 

-0.30 

C18181_0205 

SGB-GH-Weak 

0.157 

0.004 

0.204 

0.012 

5237 

3.38 

-0.51 

C18182_0207 

MS-GH-Strong 

0.051 

0.012 

0.106 

0.044 

- 

- 

- 

C18184_0442 

SGB-GH-Weak 

0.106 

0.004 

0.251 

0.010 

5202 

3.26 

-0.88 

C18187_0216 

SGB-GH-Strong 

0.169 

0.004 

0.160 

0.012 

5249 

3.41 

-0.43 

C18191_0554 

SGB-GH-Weak 

0.139 

0.004 

0.284 

0.012 

5219 

3.32 

-0.61 

C18195_0301 

SGB-GH-Weak 

0.136 

0.004 

0.232 

0.012 

5212 

3.30 

-0.65 

C18200_0241 

SGB-GH-Strong 

0.180 

0.004 

0.120 

0.014 

5232 

3.36 

-0.39 

C18200_0251 

MS-GH-Weak 

0.032 

0.010 

0.045 

0.032 

- 

- 

- 

C18200_0351 

SGB-GH-Strong 

0.189 

0.004 

0.090 

0.014 

5223 

3.33 

-0.35 

C18204_0418 

SGB-GH-Weak 

0.102 

0.004 

0.255 

0.010 

5178 

3.18 

-0.95 

C18204_0521 

SGB-GH-Weak 

0.138 

0.008 

0.185 

0.036 

5364 

3.60 

-0.50 

C18204_0524 

SGB-GH-Weak 

0.152 

0.004 

0.244 

0.012 

5215 

3.31 

-0.55 

C18206_0533 

Anomalous - Strong GN/GH 

0.117 

0.004 

0.303 

0.012 

- 

- 

- 

C18211_0559 

Anomalous - Strong GN/GH 

0.298 

0.008 

0.121 

0.044 

- 

- 

- 

C18225_0537 

Anomalous 

0.254 

0.004 

0.043 

0.012 

5212 

3.30 

-0.03 

C18225_0533 

MS-GH-Weak 

0.085 

0.014 

-0.006 

0.062 

- 

- 

- 

C18240.0320 

SGB-GH-Weak 

0.090 

0.004 

0.222 

0.010 

5190 

3.22 

-1.05 

018240.0621 

SGB-GH-Strong 

0.212 

0.004 

0.410 

0.010 

5192 

3.23 

-0.22 

018240.0705 

Anomalous 

0.198 

0.004 

0.023 

0.012 

5206 

3.28 

-0.32 

018241.0216 

SGB-GH-Weak 

0.145 

0.004 

0.272 

0.010 

5202 

3.26 

-0.60 

018243.0634 

SGB-GH-Weak 

0.087 

0.004 

0.235 

0.014 

5219 

3.32 

-1.05 

018243.0638 

MS-GH-Weak 

0.040 

0.006 

0.036 

0.020 

- 

- 

- 

018246.0716 

SGB-GH-Strong 

0.186 

0.004 

0.210 

0.012 

5206 

3.28 

-0.37 

018257.0756 

SGB-GH-Weak 

0.149 

0.004 

- 

- 

- 

- 

- 

018262.0723 

MS-GH-Weak 

0.020 

0.008 

- 

- 

- 

- 

- 

018268.0132 

SGB-GH-Weak 

0.107 

0.004 

- 

- 

- 

- 

- 

018285.0147 

SGB-GH-Strong 

0.169 

0.004 

- 

- 

- 

- 

- 

018384.0725 

SGB-GH-Strong 

0.175 

0.004 

0.357 

0.014 

5170 

3.16 

-0.44 

018384.0728 

SGB-GH-Strong 

0.217 

0.006 

0.193 

0.022 

5178 

3.18 

-0.24 









Table 3—Continued 


Star 

Comment 

I(CH) 

el(CH) 

S(3839) 

eS(3839)'^ 

Te// 

log g 

[C/Fe] 

e[C/Fe] 

C18386_0709 

SGB-CH-Strong 

0.186 

0.004 

0.303 

0.016 

5190 

3.22 

-0.37 

0.03 

C18408.0805 

SGB-CH-Weak 

0.140 

0.004 

0.340 

0.014 

5166 

3.14 

-0.64 

0.03 

018412.0233 

SGB-CH-Strong 

0.165 

0.004 

0.285 

0.012 

5226 

3.34 

-0.47 

0.02 

C18413.0647 

SGB-CH-Strong 

0.177 

0.004 

0.269 

0.020 

5245 

3.40 

-0.38 

0.02 

C18422.0256 

MS-CH-Weak 

0.021 

0.006 

0.008 

0.020 

- 

- 

- 

- 

C18422.0306 

SGB-CH-Weak 

0.111 

0.004 

0.285 

0.014 

5237 

3.38 

-0.80 

0.04 

C18422.0318 

MS-CH-Weak 

0.052 

0.006 

0.019 

0.018 

- 

- 

- 

- 

C18422.0748 

MS-CH-Weak 

0.023 

0.010 

0.041 

0.034 

- 

- 

- 

- 

C18422.0757 

SGB-CH-Strong 

0.201 

0.004 

0.052 

0.022 

5206 

3.28 

-0.31 

0.03 

C18429.0337 

SGB-CH-Strong 

0.181 

0.004 

0.263 

0.012 

5266 

3.45 

-0.34 

0.03 

C18437.0548 

MS-CH-Strong 

0.051 

0.012 

0.055 

0.040 

- 

- 

- 

- 

C18442.0658 

SGB-CH-Weak 

0.150 

0.004 

0.195 

0.016 

5239 

3.38 

-0.55 

0.02 

C18448.0557 

Anomalous 

0.218 

0.004 

0.043 

0.016 

5208 

3.28 

-0.22 

0.02 

C18461.0520 

SGB-CH-Strong 

0.189 

0.004 

0.139 

0.016 

5245 

3.40 

-0.34 

0.02 

C18478.0505 

MS-CH-Weak 

0.039 

0.006 

-0.008 

0.018 

- 

- 

- 

- 

C18478.0508 

SGB-CH-Weak 

0.115 

0.004 

0.281 

0.012 

5216 

3.31 

-0.78 

0.04 

C18484.0608 

SGB-CH-Weak 

0.149 

0.004 

0.283 

0.014 

5222 

3.33 

-0.55 

0.03 

C18488.0120 

SGB-CH-Strong 

0.212 

0.004 

- 

- 

- 

- 

- 

- 

C18496.0637 

SGB-CH-Weak 

0.108 

0.004 

- 

- 

- 

- 

- 

- 

C18501.0410 

MS-CH-Strong 

0.047 

0.008 

- 

- 

- 

- 

- 

- 

C18502.0146 

SGB-CH-Weak 

0.145 

0.004 

- 

- 

- 

- 

- 

- 

C18502.0405 

SGB-CH-Strong 

0.164 

0.004 

- 

- 

- 

- 

- 

- 

C18502.0447 

SGB-CH-Strong 

0.175 

0.004 

- 

- 

- 

- 

- 

- 


'^The uncertainty in the zero point for S(3839) is not included in the values listed above, and is ±0.019. 
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Table 4. Model Parameters, Colors, and Resulting Indices 


Te// 

log g 

Mv 

V-I 

[C/Fe] = 
I(CH) 

-0.35, [N/Fe]=+0.10 
S(3839) 

[C/Fe] = 
I(CH) 

-0.45, [N/Fe]=+1.0 
S(3839) 

[C/Fe]=-0.55, [: 
I(CH) S 

5103 

2.94 

1.887 

0.851 

0.225 

-0.001 

0.189 

0.278 

0.175 

5159 

3.12 

2.288 

0.836 

0.217 

-0.010 

0.182 

0.243 

0.167 

5216 

3.31 

2.689 

0.821 

0.208 

-0.018 

0.174 

0.209 

0.159 

5291 

3.50 

3.092 

0.802 

0.196 

-0.029 

0.164 

0.168 

0.149 

5379 

3.61 

3.293 

0.778 

0.181 

-0.042 

0.151 

0.120 

0.137 

5626 

3.79 

3.500 

0.714 

0.131 

-0.069 

0.109 

0.006 

0.098 

5997 

3.99 

3.696 

0.625 

0.070 

-0.094 

0.061 

-0.078 

0.056 

6148 

4.11 

3.889 

0.592 

0.056 

-0.102 

0.050 

-0.094 

0.047 

6209 

4.21 

4.089 

0.580 

0.052 

-0.104 

0.047 

-0.097 

0.045 

6214 

4.28 

4.287 

0.581 

0.053 

-0.103 

0.047 

-0.096 

0.045 

6186 

4.35 

4.488 

0.589 

0.055 

-0.100 

0.049 

-0.092 

0.046 

6133 

4.41 

4.648 

0.589 

0.061 

-0.095 

0.053 

-0.084 

0.050 

6062 

4.46 

4.892 

0.618 

0.070 

-0.089 

0.060 

-0.074 

0.055 

5980 

4.51 

5.094 

0.636 

0.082 

-0.083 

0.069 

-0.060 

0.063 





Table 5. 


Changes in Derived C and N Abundances for Different Model Parameters 


Star 

[C/Fe] 

[N/Fe] 

A(m- M)y 
A [C/Fe] 

= -0.10 
A [N/Fe] 

[O/Fe] 
A [C/Fe] 

= +0.15 

A [N/Fe] 

c^vc 

A [C/Fe] 

13_ 

A [N/Fe] 

[Fe/H] = - 
A [C/Fe] A 

C18200_0351 

-0.36 

0.13 



0.03 

0.01 

0.00 

-0.02 

0.08 

C18429_0337 

-0.36 

0.83 

0.03 


0.05 

0.00 

-0.08 

0.06 

0.09 

C18243.0634 

-1.08 

1.41 



0.04 

0.01 

0.01 

-0.02 

0.14 


^■Reduced from the adopted value of 10. 









